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Abstract 

First, we brieftly discuss the relationship of manifolds with density, Mean Curvature 
Flow, Ricci Flow and Optimal Transportation Theory. 

The aim of this paper is twofold. On the one hand, the study of gradient Schrodinger 
operators on manifolds. We classify the space of solutions under global conditions. As an 
application, we extend the Naber-Yau Liouville Theorem. We relate those results to weighted 
Hff,— stable hypersurfaces immersed in a manifold with density (Af, g, 4>) and we apply them 
to the classification of self-similiar solutions to the mean curvature flow. In particular, we 
give conditions for the nonexistence of self-similar solutions (shrinker or expander) to the 
mean curvature flow under global hypothesis, extending the recent results of T. Colding 
and W. Minicozzi. The Naber-Yau Liouville Theorem will be used for the classification of 
gradient Ricci solitons. We extend the Hamilton-Ivey-Perelman classification of complete 3- 
dimensional shrinking gradient Ricci solitions of bounded curvature to higher dimensions, we 
relax the hypothesis to a L 2 — type bound on its scalar curvature and an inequality between 
the scalar curvature and the Ricci curvature. 

On the other hand, the topological and geometric classification of complete weighted 
Hff,— stable surfaces immersed in a three-manifold with density (Af,g,<fi) whose Perelman 
scalar curvature, in short, P-scalar curvature, is nonnegative. Also, we classify weighted 
stable hypersurfaces in a manifold with density and a lower bound on its Bakry-Fmery-Ricci 
tensor. Moreover, we study the classification of three-manifolds with density (Af, g, <j>) under 
the existence of a certain compact weighted area-minimizing surface and a lower bound of 
its P-scalar curvature. Here, the P-scalar curvature is defined as R^ = R — 2A g <fi — |V fl (/>| 2 , 
being R the scalar curvature of (Af,g). 
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1 Introduction 

In a Riemannian manifold (N,g) there is a natural associated measure, that is, the Riemannian 
volume measure dv g = dv. More generally, we can consider Riemannian measure spaces, that 
is, triples (Af,g,m), where m is a smooth measure on N . Equivalently by the Radon-Nikodym 
Theorem we can consider triples (M,g,<p), where (j) £ C°°(M) is a smooth function so that 
dm = e^dv. The triple [M, g, (j)) is called a manifold with density (f>. 

One of the first examples of a manifold with density appeared in the realm of probablity 
and statistics, Euclidean space with the Gaussian density e -7r l x l (see [8] for a detailed expo- 
sition in the context of isoperimetric problems). In 1985, D. Bakry and M. Emery [1] studied 
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manifolds with density in the context of difussion equations. They introduced the so-called 
Bakry-Emery-Ricci tensor in the study of diffusion processes given by 

Ric^ = Ric - V 2 0, (1.1) 

where Ric is the Ricci tensor associated to (M,g) and V 2 is the Hessian with respect to the 
ambient metric g. However, manifolds with density appear in many other fields of mathematics. 



1.1 Hypersurfaces in manifolds with density 

M. Gromov [25] considered manifolds with density as mm-spaces and introduced the generalized 
mean curvature of a hypersurface E C (J\f, g, (ft) or weighted mean curvature as a natural 
generalization of the mean curvature, obtained by the first variation of the weighted area 

H <p = H + g(N,V(ft), (1.2) 

where H denotes the usual mean curvature of E and N is the unit normal vector field along E. 

Remark 1.1. In this paper, the mean curvature will be just the trace, and not the mean of the 
trace. Moreover, all the manifolds and submanifolds will be considered smooth, connected and 
orientable. Moreover, unless stated otherwise, they will be complete without boundary. 

Let (M, g, (ft) be a Riemannian manifold with density eft G C°°(M), (ft > 0. For any Borel set 
Q C JV with boundary E := 50 and inward unit normal N along E, the weighted volume of 
Q and the weighted area of Q are given by 



V*(fi):= f e*dv, A (E):= f e^dw, 
Jn Jt, 



where dv and dv% are the volume and area element with respect to g respectively. Note that 
the usual volume and area with respect to g are given when (ft = 0. 

Let ^i, t 6 (—£,£), be a smooth family of diffeomorphisms of N so that is the identity. 
Assume that 

d 



dt 



* t | E = X + uN, 



t=o 

where A is a tangential vector field along E and u is a smooth function with compact support 
on S. An immersed hypersurface E C (M,g,(ft) is weighted minimal, in short (ft— minimal, if 
it is a critical point of the weighted area functional, i.e., if 



d_ 

di 



^(^(S)) = 0, 

t=o 
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for all compactly supported variations. From Bayle's variational formulae [2] we can see that 



d_ 

dt 



MM*)) 



t=0 



that is, E is 0— minimal if and only if the weighted mean curvature vanishes. 

More generally, an immersed hypersurface E has constant weighted mean curvature = H$ 
if and only if 

(iV* t (E))-ffoW*(E))) = 

ai t=o 

for all compactly supported variations {^t}- 

Definition 1.1. We say that an open domain Q C (Af,g, <$>) with boundary E = dQ is weighted 
stationary if there exists a constant Hq so that 



d_ 

dt 



t=o 



(^(* t (E))-flo^(* t (E))) = 



for all compactly supported variations {^t}- 

As a first sight, the set of all weighted stationary domains is too big. So, we can focus on 
those which (locally) minimize the functional 

i^(* t (E)) := 4»(*t(£)) - ^(* t (E)) 

up to second order, that is, weighted stationary domains with nonnegative second variation for 
all compactly supported variations. 

Definition 1.2. We say that a weighted stationary domain Q C (M,g, 4>) with boundary E = dfl 
is stable if 

-J 0M* t (£)) - i?^(* t (s))) > o, 

/or compactly supported variations {^t}- In this case, we say that E is weighted i^— stable, 
in short, iZ^— stable. 



Here, we are adopting the strong notion of stability and not the weak version that deals 
with isoperimetric problems, i.e., a weighted stationary domain Q with boundary E := d£l is 
weighted weakly stable if 



dt 2 



t=o 



(4^(s)) - H^iMm > o, 



for all compactly supported variations {^t} that fix the weighted volume. In this case, we say 
that E is weighted isoperimetric (see [8] for details). From Bayle's variational formula [2], 
one can prove: 
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Lemma 1.1. Consider a weighted stationary open set Cl C (N,g,<j>). Let N be the inward unit 
normal along E = dCl, and H<p the constant weighted mean curvature of E with respect to N . 
Consider a compactly supported variation of CI, ^f t , with normal part uN on E. Then, we have 

j (mc%(N,N) + \A\ 2 )u 2 e' p dvj ] < J |Vu|V dvj;, (1.3) 

where \A\ is the squared norm of the second fundamental form and Ric^ 1 is the Bakry-Emery- 
Ricci curvature. 

Equivalently, we have: 

(Af - V^)"(0) = - j^uL^ue* dvx, (1.4) 

where 

L^u :=Au + g (V<f>, Vu) + (\A\ 2 + Ricf(N, N))u. (1.5) 

Here A and V are the Laplacian and Gradient operators with respect to the induced metric 
on E, A is the second fundamental form o/E, Ric is the Ricci curvature of the ambient manifold 
and V 2 is the Hessian operator with respect to the ambient metric. 

1.2 Mean Curvature Flow 

It is interesting to recall here that self-similiar solutions to the mean curvature flow in R™ +1 
can be seen as weighted minimal hypersurfaces in the Euclidean space endowed with the corre- 
sponding density (c.f. G. Huisken [18] or T. Colding and W. Minicozzi [13, 14]). We will explain 
this in more detail. 

Let X : (0, T) x S — > IR n+1 be a one parameter family of smooth hypersurfaces moving by 
its mean curvature, that is, X satisfies 




= -HN 



where TV is the unit normal along E$ = X(t, E) and H is its mean curvature. Self-similiar 
solutions to the mean curvature flow are a special class of solutions, they correspond to solutions 
that a later time slice is scaled (up or down depending if it is expander or shrinker) copy of an 
early slice. In terms of the mean curvature, E is said to be a self-similar solution if it satisfies 
the following equation 

H = ~g(x,N), (1.6) 

where c = ±1, x is the position vector in M n+1 and <?(-,•) is the standard Euclidean metric. 
Here, if c = — 1 then E is said a self-shrinker and if c = +1 then E is called self-expander. 

It is straightforward to check that self-shrinker (resp. self-expander) are weighted minimal 
hypersurfaces in (M n+1 , go, 4>-\) with density 4>-\ := — ^j- (resp. 4>+i '■= ^j - )- 
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1.3 Ricci Flow 

In G. Perelman's work [38] on the Poincare conjecture, he was able to formulate the Ricci flow 
as a gradient flow (we follow Topping's Book [43]). Let us consider a (closed) manifold with 
density (J\f, g, 4>) and introduce the following Fischer information functional 



dm, 



IN 

where dm = dv and R is the scalar curvature of (Af,g). Equivalently, integrating by parts, 



Jn 



here, is the Perelman Scalar Curvature, in short P-scalar curvature, given by 

R% = R-2A g cj>-\Vc/>\ 2 , (1.7) 

where A g and V are the Laplacian and Gradient operators with respect to the ambient metric 
g respectively. 

Remark 1.2. We focus on the compact case, i.e., M is closed. Nevertheless, the discussion can 
be extended to complete manifolds under appropiate conditions. Moreover, we will consider here 
the smooth case. 

Let us consider the variation of T which preserves the distorted volume <\>dv. So, the 
evolution of T is given by 

d F{g, 0) = - / g ( Ric - V*0, ^) eUv = - [ g (mcT, ^) e* dv. 



Therefore, if we have a solution to the coupled system 

^ - -2Ric- 
dt ~ KC ^' 



(1.8) 



(R-A 6 



then T evolves as 



^JXs^)=2^|Ric^|V(fo>0. 

Summarizing, if we define \i := dv, which is constant in time, we could view g as a gradient 
flow for the functional g — > F(g, (f), where 4> is determined by fi and dv. 
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Perelman's trick was to show that a solution of the coupled system (1.8) is somehow equiv- 
alent to a solution to the decoupled system 

( dq „ 

< f, (1.9) 

One way to see this is that solutions to the coupled system (1.8) may be generated by pulling 
back solutions of the decoupled system (1.9) by an appropiated time-dependent diffeomorphism. 
We can also reverse the argument. 

The aim of Perelman for introducing T (and the more general W— entropy) was the classi- 
fication of (gradient) Ricci solitons, that is, self-similar solutions to the Ricci flow. They can 
be analytically described as manifolds with density (J\f, g, <j>) such that Ric^ 1 = Xg, with Ael. 
They are called shrinking, steady or expanding depending if A > 0, A = or A < respec- 
tivelty. Steady Ricci solitions appear as critical points of the T funtional (shrinking Ricci solitons 
appear as critical points of the W— entropy functional). See [43] for a detailed exposition. 

Remark 1.3. Perelman •proved that is not the trace of Ric^ 1 but they are related by a 
Bianchi type identity. Moreover, Perelman's entropy functional and manifolds with density have 
applications to String Theory (see [6]). 



1.4 Optimal Transportation Theory 

Recently, in Optimal Transport Theory (we follow Part II in Villani's Book [44]), geometric 
problems on manifolds with density (jV, g, <p) has been related to correspondent problems in the 
Wasserstein space of probability measures equipped with the quadratic Wessertein metric, which 
corresponds to a relaxed version of Monge's optimal transportation problem. We will explain 
the relationship between them from a formal point of view, we will not be rigourous (see [44] 
and references therein). Moreover, we will assume all the objects involved are smooth and verify 
the appropiate convergence conditions to avoid technical problems and clarify the exposition. 

Let (M,g) be a Riemannian manifold. Denote by dv and d the Riemannian meassure and 
distance with respect to g respectively. The Wasserstein space T > 2{-N') is the space of probability 
meassures which have a finite moment of order 2, i.e., 

V 2 {M) = j/x G V{M) : j d( Po ,p) 2 dfi < +oo 

where po G N is arbritrary and V{N) is the space of probability meassures. The Wasserstein 
space will be equiped with the Wasserstein distance 

W 2 (/x, v) = ^inf | J d(p, q) 2 dir : n G 1%, v) 
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where the infimum is taken over all joint probability meassures ir on M x M with marginals p 
and v. Such joint meassures are called transference plans; those achieving the infimum are 
called optimal transference plans. 

Remark 1.4. We should think on d 2 as the cost function on the optimal transportation problem. 

F. Otto [35] noticed that computations of Riemannian nature can shed light on Optimal 
Transportation Theory. So, the problem was to establish rules for formally perform differential 
calculus on ^(AT)- We focus here on a certain class of functional that appear naturally on the 
theory, as for example, the Boltzman entropy. 

Let us consider <f> : N — > IR a smooth function to distort the reference volume meassure, 
i.e., dm = e^ dv, that is, we can consider a manifold with density (M,g,(f>). Consider a smooth 
funtion U : IR + — > IR which will relate the values of the density of the probability meassure and 
the value of the functional, i.e., 

U m {p) = I U(p(p))dm, p = pdm. 

Remark 1.5. We should think on U m as the internal energy of a fluid, i.e., it is like the energy 
contained on a fluid of density p. The function U should be thought as a property of the fluid 
and it might mean some microscopic interaction. It is natural to assume U(0) = 0. 

In analogy with thermodynamics, we can introduce the preassure as: 

p(p) :=pU'(p)-U(p). 

Remark 1.6. The physical meaning of the preassure says that if a fluid is enclosed in a domain 
Q, then the preassure felt by the boundary dfl at a point is normal and proportional to p at that 
point. 

So, one can consider the total preassure 




and again, compute the variation of this functional with respect to small variations of the 
meassure, which leads to the iterated preassure: 

P2(p) ■= PP'(p) - P(p)- 

One can see that the preassure and iterated preassure appear when we differentiate the energy 
functional U m to first and second order respectively. F. Otto [35] gave an explicit expression for 
the gradient and Hessian of the functional U m . For a given meassure p, the gradient of U m at p 
is a tangent vector at p in the Wasserstein space, and it is given by 

grad U m = - (Ap(p) + g(V(/),Vp(p))) dm, 
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where, remember, dm = dv. The Hessian of U m at p, is a quadratic form on the tangent space 
TuV2{N), but its expression is rather complicated in general. We continue this discussion only 
for the Boltzman entropy funtional, which is given by 

M-min) '■= / pin pdm, p = pdm. 

For the Boltzman entropy, its gradient is given by 

grad M % m = - (Ap + g(V<p, Vp)) dm 
so, we can see that critical points of the Boltzman entropy correspond to positive solution to 

Ap + g(V<j),Vp) = 0, 

or, equivalently, they correspond to positive Jacobi functions for the gradient Schrodinger oper- 
ator 

L^p:= Ap + g(V0,Vp). 

2 Organization of the paper and Main results 

We will describe here the organization of the paper and we will outline the main results. More- 
over, we should remark that the required differentiability is much lower of that we consider here, 
we prefer to avoid technical details assuming smoothness for clarifying the exposition and ideas 
involved in this paper. 

In Section 3 we develop the theory of gradient Schrodinger operators defined on a manifold 
S acting on piecewise smooth functions of compact support u G (7q°(S). That is, given (E,g) 
a Riemannian manifold, complete or compact with boundary possibly empty, (p G C°°(S) and 
q € C°°(£) called potential, we study differential operators on the form 

L^u := An + g(V(f>, Vu)+qu, u G Cg°(E). (2.1) 

When V</> = 0, that is, we are dealing with an usual Schrodinger operator L := A + q, it 
is known that L is self-adjoint with respect to the L 2 — inner product, but this is not longer 
true for a gradient Schrodinger operator. Nevertheless, we will see that is self-adjoint with 
respect to a weighted L 2 — inner product, the weight depends on <j>. We will show that gradient 
Schrodinger operators mimic most of the properties of the usual Schrodinger operators when we 
consider this weighted L 2 — inner product. When S is compact and there exists a non identically 
zero solution to L^u = 0, then u vanishes nowhere and the linear space of such functions is one 
dimensional. We would like to extend this property to the noncompact case impossing a global 
condition on the triple (X,<7, 4>). 
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Definition 3.2. Let (£,g) be a complete manifold and <\> G C°°(S). We say i/iai 
(E,g,4>) is of finite type if there exists a sequence of cut-off functions C 
such that 

• < ipi < 1 m S. 

• The compact sets 0^ := form an increasing exhaustion o/S. 

• XTie sequence of weighted energies {J s |V^i| 2 e^ dfEj^ is bounded. 

Moreover, given a complete manifold (£,<?) a^c? smooth functions <fi,q G C 00 (S) ; we 
say i/iai i/ie gradient Schrddinger operator with potential q given by (2.1) zs o/ 
finite type if(T,,g,(p) is of finite type. 

In the case V</> = 0, Definition 3.2 says that (£, g) has finite capacity (see [24]). The following 
result, which extends [28, Theorem 2.3], tells us that finite type gradient Schrodinger operators 
behaves as in the compact case. 

Theorem 3.2. Let (£,<?) be a complete manifold and take (f>,qi,q2 G C°°(£). As- 
sume that 

1. The gradient Schrddinger operator £<£ )?1 := A + g(V/, V-) + q\ is of finite type 
and there is a positive subsolution u, i.e., u satisfies 

L<p >qi u < 0. 

2. There exists a bounded function v G C°°(E) so that 

L^v := Av + g(Vf, Vv) + q 2 v > 0. 

3. qi - 02 > on S. 
Then, v/u is constant. 

Moreover, we above technique allows us to extend the Naber-Yau Liouville Theorem [33]: 

Theorem 3.3. Let (£,g) be a complete manifold and take 4>, q G C°°(S). Assume 
the gradient Schrddinger operator L^ q := A + g(V(fr, V-) + q is of finite type. 

1. Assume q > 0. // there is a nonnegative (nonidentically zero) subsolution u, 
i.e., u satisfies L$ A u < 0, then u is constant and either u = or q = on S. In 
particular, any solution bounded above or below to Lf q u = must be constant. 

2. Assume q < and L^ A is stable. If there is a bounded supersolution u, i.e., u 
satisfies L<p >q u > 0, then u is constant and either u = or q = on S. 
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3. Assume q < and there is a nonnegative supersolution u, i.e., u satisfies 
L<h qU > 0, such that the sequence is uniformly bounded. 

Then u is constant and either u = 0orq = 0onT,. Moreover, if q = 0, we 
only need to assume that u is bounded below. 

In Section 4, we will apply Theorem 3.2 to recover and extend a recent result of Colding- 
Minicozzi [13, 14] on the nonexistence of certain self-shrinkers. To understand properly the main 
results in this section, we need the following definition: 

Definition 4.2. Let S C IR ra+1 be a self-shrinker (resp. expander). We say that £ 
is of finite type if (£, g, 4>-\) (resp. (T,,g,(f) + i)) is of finite type. Here, g denotes 
the induced metric on S. 

Here, <j) c = c^j-, where x is the position vector in IR n+1 , and c = — 1 if £ is a self-shrinker 
and c = +1 if S is a self-expander. T. Colding and W. Minicozzi [13, 14] proved that there 
are no L— stable self-shrinker £ with polynomial area growth. We say that a self-shrinker (resp. 
expander) £ is L— stable if it is stable as a weighted minimal hypersurface in (M ra+1 , go, <fi-\) 
(resp. (R n+1 , go, <fi+i)). Polynomial area growth implies that 

f e^ 1 dv -»• 0, 

JB(i+l)\B(i) 

where B{r) denotes the geodesic ball of radius r centered at a fixed point. In particular, S is of 
finite type. So, we can prove: 

Theorem 4.1. There are no complete L— stable self-shrinkers in M n+1 of finite type. 

As far as we know, there are no general classification results for L— stable self-expanders. 
We can, at least, give some condition for the nonexistence: 

Proposition 4.1. There are no complete L— stable self- expanders in M n+1 of finite 
type and \ A\ 2 > 1/2. 

In fact, there is another notion of stability more related to graphs. Let S C IR n+1 be a 
self-shrinker (resp. expander). Let a G M n+1 be a constant vector and consider u a := g(N,a), 
then u a satisfies (see [13, 14]): 

L u a := Au a + g (V0 C , Vu a ) + \A\ 2 u a = 0. 

If S is a complete multi-graph, say with respect to the e n+ \ direction, then v := g (iV, e n+ \) 
is positive and satisfies the above equation, therefore Lemma 3.1 yields that Lo is stable in the 
sense of Definition 3.1. Therefore, a self-similar solution to the mean curvature flow is Lo— stable 
when Lo defined above is stable in the sense of Definition 3.1. Thus, this leads us to: 
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Theorem 4.2. The only complete self-similar solutions to the mean curvature flow 
that are Lq— stable and of finite type are hyperplanes. 

In particular, The only self-similar solutions (shrinker or expander) to the mean cur- 
vature flow that are complete multi-graphs of finite type are hyperplanes. Moreover, 
the only self-shrinkers that are entire graphs are the hyperplanes. 

Recently, L. Wang [45] proved that the only self-shrinkers that are entire graphs are hyper- 
planes. This result follows from us since a proper self-shrinker has polynomial volume growth 
[11, Theorem 1.3], and hence, it is of finite type. 

In Section 5, we use Theorem 3.3 for classifying gradient Ricci solitons. Noncompact gradient 
Ricci solitons are fundamaental in the proof of the Poincare Conjecture. From the works of R. 
Hamilton [26], T. Ivey [27] and G. Perelman [38, 39] we have: The only three dimensional 
shrinking gradient Ricci solitons with bounded curvature are the finite quotients q/IR 3 ; S 2 x H. 
or S 3 . 

P. Petersen and W. Wylie [36] showed that the only complete shrinking gradient Ricci soliton 
such that J* E [Ricpe^ dv-z < oo and vanishing Weyl Tensor, i.e. W = 0, are finite quotients of 
JBL n , S n_1 x ffi or S" (see [36, 37] and references therein for a detailed exposition). 

We extend the above result imposing only a L 2 — type bound on the scalar curvature R. 
Moreover, they assume that the Weyl tensor vanishes identically, we change that condition by 
an inequality between the scalar curvature and the Ricci curvature. 

Theorem 5.1. Let (S,<7, <f>, A) be a complete shrinking soliton of dimension n > 3 
so that XR < |Ric| . A ssume that there exists a sequence of cut-off functions 
{V'jjj C Cq°(E) such that 0<^<lmE, the compact sets VLi := ■0~ 1 (1) form an in- 
creasing exhaustion ofH, and the sequence of weighted energies {/ s \Vi^i\ 2 R 2 e^ dvj^^ 
is bounded. 

Then, (E,g) has constant scalar curvature and XR = |Ric| 2 . 

So we recover the Hamilton-Ivey- Perelman classification assuming only a L 2 — type bound on 
the scalar curvature and the inequality between the scalar curvature and the Ricci tensor. We 
also give the following classification result when the scalar curvature does not change sign: 

Theorem 5.2. 

• A complete gradient expander soliton of finite type and nonnegative scalar cur- 
vature is Ricci flat. 

• A complete gradient steady solition of finite type whose scalar curvature does 
not change sign is Ricci flat. Moreover, if <f> is not constant then it is a product 
of a Ricci flat manifold with R. Also, £ is diffeomorphic to W 1 . 
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• A complete gradient shrinking soliton of nonpositive scalar curvature is Ricci 
flat. 

In Section 6 we introduce the key tool for studying i?^— stable surfaces. We will associate a 
self-adjoint Schrodinger operator to any weighted H^— stable surface which is stable in the sense 
of Schrodinger operators (see Definition 6.1), i.e., its first eigenvalue is nonnegative. Specifically, 
we prove 

Lemma 6.1. Let S C (A^ 3 , <7, 0) be a complete weighted H^— stable surface. Then, 

[ {V - aK)f 2 dw < [ \Vf\ 2 dvx 

JT, JY, 

for any f G Cq°(£) compactly supported piecewise smooth function, where 

In other words, the Schrodinger operator 

L:=A-aK + V 

is stable in the sense of Definition 6.1. 

Section 7 is devoted to the classification of complete weighted H$— stables surfaces in a 
manifold with density <p € C°°(S) so that its P-scalar curvature is nonnegative. First, we 
extend the Rosenberg's diameter estimate for weighted -f/^,— stable surfaces given a positive 
lower bound for + (see [40] for the Riemmanian case), 

Theorem 7.1. Let X C (N,g,4>) be a weighted H^— stable surface with boundary 
as. // K» + Hj> c> on E, then 

2tt 

ds (p, d'E) < —= for all p £ S, 
V3c 

where d^ denotes the intrinsic distance in E. Moreover, if S is complete without 
boundary, then it must be topologically a sphere. 

Theorem 7.1 says that the only complete noncompact weighted H^— stable we shall consider 
are the minimal ones. So, the next step is to classify the complete noncompact weighted minimal 
stable surfaces in the spirit of the works of D. Fischer- Colbrie and R. Schoen [23] and R. Schoen 
and ST. Yau [41]. The theorem we provide here extends a previous result of P. T. Ho [17]. In 
[17], the author needs the additional hypothesis that |V0| is bounded on S. Here, we drop that 
condition and characterize the stable cylinder as totally geodesic, flat and with constant density 
along S. We show: 
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Theorem 7.2. Let E C (N,g,(ft) be a complete (noncompact) weighted stable mini- 
mal surface where > 0. Then, E is conformally equivalent either to the complex 

plane C or to the cylinder § 2 xl. In the latter case, E is totally geodesic, flat and 
(ft is constant and = along E. 

We finish Section 7 classifying ify— stable hypersurfaces in a manifold with density satisfying 
a lower bound on the Bakry-Emery-Ricci tensor. Specifically: 

Theorem 7.3. Let E C (Af,g,4>) be a complete H^— stable hypersurface of finite type 
and assume that Ric^ 1 > k, k > 0. Then, E is totally geodesic and Ric^(iV, N) = 
a/ong E. 

In particular, Theorem 7.3 extends results obtained by X. Cheng, T. Mejia and D. Zhou [10]. 

In Section 8 we give local splitting results for three-manifolds with density (M, g, (ft) verifying 
some lower bound on its P-scalar curvature R^ . The idea is to extend the splitting theorems 
developed by Cai-Galloway [5], Bray-Brendle- Neves [4] and I. Nunes [34]. 

Theorem 8.1. Let (M 3 ,g,(ft) be a complete manifold with density containing a 
compact, embedded weighted area-minimizing surface E. 

1. Suppose that Roo > Xe^, for some positive constant X, and A^iY,) = Then 
E has genus zero and it has a neighborhood in M which is isometric to the 
product § 2 x (— e, e) with the product metric g + \ +dt 2 (up to scaling the metric 
g); here g + \ is the metric of constant Guassian curvature +1. Moreover, (ft is 
constant in U. 

2. Suppose that R^ > and E has genus one. Then E has a neighborhood in M 
which is flat and isometric to the product T 2 x (— e, e) with the product metric 
go + dt 2 ; here go is the metric of constant Guassian curvature 0. Moreover, (ft 
is constant in U. 

3. Suppose that Roo > —Xe^, for some positive constant X, E has genus 7 > 2 and 
^(E) = 47r (^~ 1 ) _ Then E has a neighborhood in N which is isometric to the 
product E x (— e,e) with the product metric g^\ + dt 2 (up to scaling the metric 
g); here g-\ is the metric of constant Guassian curvature —1. Moreover, (ft is 
constant in U. 

Finally, in Section 9 we establish the classification results for manifolds with density (ft 
satisfying a lower bound on its P-scalar curvature and the existence of certain compact weighted 
area-minimizing surface in its homotopy class. First, we focus in the case the P-scalar curvature 
has a positive lower bound. We prove: 
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Theorem 9.1. Let (J\f 3 ,g,cf)) be a compact manifold with density so that is 
positive and tt2(M) / 0. Define 

A{M, g, 4>) :=inf{A (/(§ 2 )) : / G f) , 

where T is the set of all smooth maps f : S> 2 — > N which represent nontrivial elements 
in 7T2 (M) ■ 

Then, we have 

( R°°(x) ) 
A{M, g, 0)min j : x € M \< 8tt. 

Moreover, if equality holds, the universal cover of (Af,g) is isometric (up to scaling) 
to the standard cylinder S 2 x H. and <j) is constant in M. 

Second, when the P-scalar curvature is nonnegative: 

Theorem 9.2. Let (AA 3 , g, 0) be a complete manifold with density so that is 
nonnegative. If (J\f, g, (f>) contains a weighted area miminizing compact surface in its 
homotopy class of genus greater than or equal to 1, then the product manifold T 2 xl, 
where T 2 is a torus equiped with the standard flat metric, is an isometric covering 
of{M,g), and 4> is constant in N . In particular, (M,g) is flat. 

And finally, when the P-scalar curvature has a negative lower bound: 

Theorem 9.3. Let (M 3 ,g,cf)) be a complete manifold with density so that > 
— Ae^ for some positive constant X. Moreover, suppose that £ C (Af,g,<p) is a two- 
sided compact embedded Riemannian surface of genus 7 > 2 which minimizes area 
in its homotopy class. Then, 

Am > 4 -^. 

Moreover, if equality holds, the product manifold S 7 x R, where S 7 is a compact 
surface of genus 7 equiped with a metric of constant Gaussian curvature —1, is an 
isometric covering (up to scaling) of (Af,g), and (f> is constant in N . 

Remark 2.1. Most of the above results can be generalized to the case that S has finite index, 
that is, there is only a finite dimensional space of normal variations which strictly decrease 
the weighted area (similar to results of Fischer- Colbrie [22]). Moreover, we could relax the 
hypothesis on the nonnegativity of the P-scalar curvature by either an integrability conditon or 
a decay condition as in [3] or [20]. 
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3 Gradient Schrodinger operators of finite type 

We will follow the references [9, 15, 16]. Consider (T,,g) a compact Riemannian manifold with 
boundary <9£ (possibly empty). Set q G C°°(£) and X G X(S) a smooth vector field along E. 
Consider the differential linear operator, called generalized Schrodinger operator, given by 

L: Cg°(S) -»• C°°(£) 

u — > Lu := Au + g(X,Vu) + qu, 

where A and V are the Laplacian and Gradient with respect to the Riemannian metric g 
(respectively) and Cg°(£) stands for the linear space of compactly supported piecewise smooth 
functions on S. We also denote by |-| the norm with respect to the Riemannian metric g. 
Moreover, 

• In the case that X = V</> for some smooth function <j), i.e., L^u := Au + g(V</>, Vu) + qu, 
we call Lfj, by gradient Schrodinger operator. 

• In the case that X = 0, i.e., Lu := Au + qu, we call L simply by Schrodinger operator. 

Note that the generalized Schrodinger operator L^, given by (2.1), is a gradient Schrodinger 
operator for X := V0. In general, a gradient Schrodinger operator is not self-adjoint with respect 
to the L 2 — inner product because of the first order term, but it is self-adjoint with respect to a 
weighted inner product. 

Proposition 3.1. Let (£,#) be a compact Riemannian manifold with boundary dT, (possibly 
empty). Set q,<fi G C°°(S) and consider 

L^u := Au + g(V4>, Vn) + qu , u£ C£°(£). 

Then, L^ is self-adjoint acting on u G Co°(S) with respect to the weighted inner product 
given by 

g{u,v)^ := J uve^dvs. 
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Proof. We have to prove that 

g^L^v)^ =g(L (P u,v) (t) , for all n, v G C£°(£). 
Take u,v G C£°(£), then 

<7 (n, L^v)^ = I u L^v e^ dv-£ = I [ue^ Av + ue^giVcj), Vv) + quve^j dv^ = 

= y ^n div^Vv) + quve^j dv^ = j (—g(Vu,Vv) + guv) dvs = 

where we have used the Divergence Theorem on the second line. □ 

Let be a gradient Schrodinger operator. We say that A G 1 is an eigenvalue of if 
there is a not identically zero function n 6 C^°(E) so that 



In the case of gradient Schrodinger operators, standard spectral theory (see [15, Pages 335- 



Theorem 3.1. Let (X,g) be a compact Riemannian manifold with boundary <9E (possibly empty). 
Set q, 4> G C°°(S) and consider 



Then, 

• L^ has real eigenvalues \\ < A2 < . . . with — > +00 as k — > +00. 

• There is an orthonormal basis {n^} C Cq°(E) for the weighted L 2 metric given in Propo- 
sition 3.1 so that L^Uk = —XkUk- 

• The lowest eigenvalue \\ is characterize by 



• Any eig enf unction for \\ does not change sign and, consequently, if n G (7^(2) is another 
solution to L^u = —\\u, then u is a constant multiple of u\. 



Lcj>u = —An. 



336]) gives: 



L^u := An + g(V(f>, Vn) + qu , u G C£°(£). 
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In the case that (X,g) is complete noncompact, (j) £ C°°(S), there may not be a lowest 
eigenvalue for L^, however, we can still define the bottom of the spectrum (we still call it Ai) 
by 

f J s (|Vn| 2 -quAeUvx ) 
W>*) = ™{ V ku ^dL = « e CS-(E)\ { 0}j, 

where now the infimum is taken over piecewise smooth functions of compact support and we 
must allow Ai(£,</>) = — oo. Another way to calcule Ai (£,(/>) is 

Ai(S,(/>) = inf{Ai(^,0) : SlcE), 

where the infimum is taken over all relatively compact domains C S with (at least) C 1 
boundary. 

Definition 3.1. Let (£,<?) fre o complete manifold, (j) 6 C°°(E) and a gradient Schrddinger 
operator acting on u £ Cq°(S). We say £/ta£ is stable z/Ai(S,<^) > 0. 

We continue characterizing stable gradient Schrodinger operators by a variation of an argu- 
ment of Fischer-Colbrie and Schoen. T. Colding and W. Minicozzi [13, 14] proved a particular 
case of the following result in the context of stable self-shrinkers. We will prove here the general 
version. 

Lemma 3.1. Let (T*,g) be a complete manifold and (j) £ C°°(S). Then, the following statements 
are equivalents: 

(1) is stable. 

(2) There exists a positive function u on £ so that L^u = — Ai(£,<^)u ; Ai > 0. 

Proof. (1) implies (2): We argue as in [13, Lemma 9.25]. Since L$ is stable, Definition 3.1 
yields that Ai (£,</>) > 0. 

Fix a point p 6 S and let B(rk) be the geodesic ball in £ centered at p of radius r k . We can 
consider an increasing sequence {r k } so that dB(r k ) is at least C 1 and r k — >■ +oo as k — > +oo. 

It is clear that \\{B{rk),4>) > Ai (£,</>) > an d 

lim AiOB(r fc ),0) = Ai(E,0)>O. 

From Theorem 3.1, there exists a positive Dirichlet function u k on B(r k ) so that 

L<t,u k = -\i(B{r k ),<j>)u k , 
so that (after multiplying by a constant) u k (p) = 1. This holds for all fc. 
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Then, the sequence {uk} is uniformly bounded on compact sets of £ by the Harnack In- 
equality [16, Theorem 8.20]. Also, {uk} has all its derivatives uniformly bounded on compact 
subsets of S by Schauder estimates [16, Theorem 6.2]. Therefore, Arzela-Ascoli's Theorem and 
a diagonal argument give us that a subsequence (that we still denote by {uk}) converges on 
compact subsets of S to a function u G C°°(E) which satisfies 

' L^u = — Ai(S, 4>)u in S 
u > in £ 

= 1 at p£ S, 

again, by the Harnack Inequality, we obtain that u > in £. This proves (1) implies (2). 

(2) implies (1): We argue here as in [14, Proposition 3.2]. Let u be the positive smooth 
function given by item (2). Then, it satisfies L<pu < since Ai(£, <j>) > 0. Set w = lnu, then 

Aw = | Vw | 2 < —g(V(f>, Vw) — q — \Vw\ 2 , 

which implies 

div(e^W) < - (q+ \w\ 2 ) e*. 



Given tp G C^°(£), applying Stokes' Theorem to div (if> 2 e^V w) gives 
= j (^e^V^, Vw) + V 2 div ^Vw) ) cfo s 

< J (2^g(V^,Vw)-^ 2 (q+\Vw\))e' t> dvj: 

< J (iVVl 2 -^ 2 ^) e*dvz 

where we have used 2ipg(Vip, Vw) < \Vip\ 2 + V 2 \Vw\ 2 . Therefore, Ai(£, <j>) > 0. □ 

We can see from Theorem 3.1 that, if £ is compact and there exists a non identically zero 
solution of L^u = 0, <j> G C°°(£), then it vanishes nowhere and the linear space of such functions 
is one dimensional. We would like to extend this property to the noncompact case impossing a 
global condition of the triple (£,#, </>). 

Definition 3.2. Let (£,g) be a complete manifold and 4> G C°°(£). We say that (E,g,(j)) is of 
finite type if there exists a sequence of cut-off functions {ipi},i C Co°(£) such that 

• < tpi < 1 in £. 

• T/ie compact sets 0^ := ~ 1 ( 1 ) /orm an increasing exhaustion o/£. 
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• The sequence of weighted energies {J E \ Vipi\ 2 e^ dv^^ is bounded. 

Moreover, given a complete manifold (£,<?) and smooth functions (f),q € C°°(S) ; we say i/tai 
£/ie gradient Schrddinger operator with potential q, given by 

L^u := Au + g(V(p, Vu) + qu , ii £ C£°(E), 

is 0/ finite type if (Y,,g,<f>) is of finite type. 

We should point out that the above definition does not depend on the potential q € C°°(S). 

Remark 3.1. In i/ie case V(f> = 0, i.e, when L is simply a Schrddinger operator, Definition 
3.2 says that (S, g) has finite capacity (see [24])- Moreover, we should remark that the required 
differentiability is much lower of that we consider here. 

The following result, which extends [28, Theorem 2.3], tells us that finite type gradient 
Schrodinger operators behaves as in the compact case. 

Theorem 3.2. Let (£,<?) be a complete manifold and take 4>,qi,q 2 £ C°°(S). Assume that: 

1. The gradient Schrddinger operator L^,^ qi := A + g(V(j), V-) + gi is of finite type and there 
is a positive subsolution u, i.e., u satisfies 

L^ qi u < 0. 

2. There exists a bounded function v £ C°°(£) so that 

L^ q2 v := Av + g(V<j>, Vu) + q 2 v > 0. 

3. qi - q 2 > on S. 
Then, v/u is constant. 

Proof. Let u,v G C°°(£) as in the statement of the theorem. A straightforward computation 
shows that 

2 

-div (u 2 e*V(v/uj) = vdivte^Vv) - — div(e^Vu). (3.1) 
u \ J u 

Now, since u and v satisfies the diferential inequalities in the statement of the theorem, 

multiplying by e^, we get 

div(e^Vu) < -q^ue^ and div(e*Vu) > -q 2 ve (f> , 

with this information and (3.1), we obtain 

^div (uVV(v/u)) > (gi - q 2 )v 2 e <t> > 0. (3.2) 
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Consider the sequence {V'ili C Co°(E) given by Definition 3.2. Then, 
= y div ((Y>f ^)(tiVV(v/u))) dvj: 

>2 / V*«we^(V^i,V(u/«))d«E+ / Vi« 2 e*|V(u/ti)| 2 cfo E , 
Je is 

that is, using Holder Inequality, 

/ ^?uV |V(w/n)| 2 dv s < -2 / ipiuve^g (V^, V(u/u)) efo E 

= -2 / ^ittve*^ (V^i, V(«/u)) 

J supp(ipi)\Qi 

< 2 ( / ^?uV |V(t;/u)| 2 cfo E J ( / vV |V^| 2 dv E 



<2[ I ^fu 2 e*\V{v/u)\ z dv^j (J vV|V^r ) • 



- i 

I WW,.. I 2 



Now, since L^ qi is of finite type and v is bounded, there exists a constant C so that 

y vV (VV'il 2 d«E < C for all i, 
which implies from the previous inequality that 

y \V(v/u)\ 2 cfo E < 4C for all i, 

therefore, we obtain 

y uV |V(?;/n)| 2 dvx < 4C (3.3) 

lim / ^ 2 u 2 e^\V(v/u)\ 2 dv E = 0. (3.4) 

Thus, (3.3) and (3.4) imply that 

/ u 2 e^\V{v/u)\ 2 cfc E < 2VC ( / V 4 Ve^|V(7;/n)| 2 dv s ) 

as i goes to +oo, thus v/u is constant on S. □ 



and 
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As a consequence of the above technique, we can extend the Naber-Yau Liouville Theorem 
(see [36]). 

Theorem 3.3. Let be a complete manifold and take (f>,q € C°°(£). Assume the gradient 

Schrddinger operator L^^ := A + g(V(f>, V-) + q is of finite type. 

1. Assume q > 0. If there is a nonnegative (nonidentically zero) subsolution u, i.e., u satisfies 
L(j>,qU < 0, then u is constant and either u = or q = on S. in particular, any solution 
bounded above or below to Lf jQ u = must be constant. 

2. Assume q < and L^ i? is stable. If there is a bounded supersolution u, i.e., u satisfies 
L<f>,qU > 0, then u is constant and either u = or g = on £. 

5. Assume q < and £/iere is a nonnegative supersolution u, i.e., u satisfies L^^u > 0, such 
that the sequence is uniformly bounded. Then u is constant and 

either u = or q = on E. Moreover, if q = 0, we on/y need to assume that u is bounded 
below. 

Proof. Let us prove each case: 

1. Assume q > 0. A straightforward computation, as above, shows that 

= J div (V-e*- u Vn) <to E 

= y ^ 2 e~"div ^V«)* s + y s(v($e- u ),e*Vu) <to E 

<2 f iJie*- u g(yiPi,Vu)dvz- f tpfe^ |Vn| 2 dt> s - / ifie^uqdw, 

therefore, using Holder Inequality and u > 0, we obtain 

ijfe^- u \Vu\ 2 dvs < [ ^e^- u \Vu\ 2 dv s + [ ^fe^uqdv^ 

i i 



< 2 / itf e*-u | Vu |^ ^ i / e ^-« | V ^.|^ 

\J supp(ipi)\Vli I \J supp(tpi)\Cli 



l l 

iff. \ 2 



< 2 ( / $e*~ w |Vn| 2 cfc s ] ( / |V^| 2 d« s 

Thus, arguying as Theorem 3.2, we obtain that u must be constant and either u = or 
9 = 0. 
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2. Assume q < and u is bounded. Then, Theorem 3.2 implies that u > 0. Therefore, 
calculating as above: 

= J div (ipfue^uj dv-s 

= J ipfu div (e^Vuj dv^ + J e^g (V(^ 2 u), Vu) cfo s 

>2 / ip i ue' t> g{Vipi,Vu)dv s + / ^ e <j> | Vu| 2 dv^ - \ ipfe^u 2 qdvs, 

that is, 

I ^e^\Vu\ 2 dw - f ifieWqdvE < -2 f i/>iue*g (VV»i,Vu) dus, 

and we can finish as in the previous case. 

3. We argue as in item 2. That is, we have 

f ^V|Vii| 2 d«E- / ifieWqdvz < -2 f faue+g (VV»i,V«) cfo s , 
which implies 

/ | Vu| 2 (fo s < 2 ( / ^?e* | Vu| 2 dv E J If v?e* |V^| 2 <fo S J , 

therefore, we can finish as in Theorem 3.2. 

□ 

Remark 3.2. Theorem 3.2 and Theorem 3.3 do not follow from [28] by doing a conformal 
change of metric. 

4 Applications to Mean Curvature Flow 

We will apply Theorem 3.2 to recover and extend a recent result of Colding-Minicozzi [13, 14] 
on the nonexistence of certain self-shr inkers. Set c = ±1, as we pointed out, a self-similar 
solution £ to the mean curvature flow satisfies H = — |(x, N), when c = — 1, S is a self-shrinker 
and when c = +1, £ is a self-expander. A self-shrinker (resp. expander) £ C IR n+1 can be 
considered as a weighted minimal hypersurface in the manifold with density (M n+1 , (•,•), (f>-i) 
(resp. (IR n+1 , (■,-), (p+i)), where (-,-) denotes the standard Euclidean metric and (fr-i := — ^j- 
(resp. 0+i := ^j-). So, 
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Definition 4.1. We say that a self-shrinker (resp. expander) E is L— stable if it is stable as a 
weighted minimal hypersuface in (R n+1 , (•, •), <p-i) (resp. (IR n+1 , (-, •), <fi+i) )■ 

We also need the following definition: 

Definition 4.2. Let E C IR n+1 be a self-shrinker (resp. expander). We say that E is of finite 

type if (E, g, <p-i) (resp. (T,,g,<fi + i)) is of finite type. Here, g denotes the induced metric on E. 

So, we can prove: 

Theorem 4.1. There are no complete L— stable self-shrinkers in WL n+1 of finite type. 

Proof. Let E C (R n+1 , (•,■), _i) be a self-shrinker solution, i.e., it satifies H = ^(x,N). More- 
over, we are assuming E is L— stable, that is, it is stable as a weighted minimal hypersurface, 
therefore the gradient Schrodinger operator (see Lemma 1.1) given by 

L:= A + g(V^ 1 ,V-) + (\A\ 2 + ±), 

is stable in the sense of Definition 3.1. From Lemma 3.1, there exists a positive subsolution, 
that is, there exists u > so that 

Lu < 0. 

Moreover, q := \A\ 2 + \ > \ > 0. Thus, Theorem 3.3 says that u is constant and q = 0, 
which is a contradiction. □ 

T. Colding and W. Minicozzi [13, 14] proved that there are no L— stable self-shrinker E with 
polynomial area growth. Polynomial area growth implies that 

/ 0_i dv ->■ 0, 

JB(i+l)\B(i) 

where B(r) denotes the geodesic ball of radius r centered at a fixed point. In particular, E is of 
finite type. 

As far as we know, there are no general classification results for L— stable self-expanders. 
We can, at least, give some condition for the nonexistence: 

Proposition 4.1. There are no complete L— stable self- expanders in M n+1 of finite type and 
\A\ 2 > 1/2. 

In fact, there is another notion of stability more related to graphs. Let E C W l+1 be a 
self-shrinker (expander). Let a £ M n+1 be a constant vector and consider u a := (N,a), then u a 
satisfies (see [13, 14]): 

L u a := Au a + g (V0 C , Vu a ) + \A\ 2 u a = 0. (4.1) 

If E is a complete multi-graph, say with respect to the e n+ \ direction, then v := (N, e n+ \) is 
positive and satisfies (4.1), therefore Lemma 3.1 yields that Lq is stable in the sense of Definition 
3.1. This motivates the following: 
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Definition 4.3. We say that a self-similiar solution to the mean curvature flow S is Lq— stable 

if the gradient Schrddinger operator 



is stable in the sense of Definition 3.1. 

Therefore, the above discussion leads us to: 

Theorem 4.2. The only complete self-similar solutions to the mean curvature flow that are 
Lq— stable and of finite type are hyperplanes. 

In particular, The only self-similar solutions (shrinker or expander) to the mean curvature 
flow that are complete multi-graphs of finite type are hyperplanes. Moreover, the only self- 
shrinkers that are entire graphs are the hyperplanes. 

Remark 4.1. It would be interesting to know is a self-shrinker that is a complete multi-graph 
must be an entire graph. 

5 Applications to gradient Ricci solitons 

As we said at the Introduction, gradient Ricci solitions (steady and shrinking) appear as critical 
points of the functional T (or W— entropy) introduced by G. Perelman [38]. Gradient Ricci 
solitions are self-similar solutions of the Ricci flow and they arise as limits of dilatations of 
singularities of the Ricci flow. Analitically, we can define a Ricci solition as a manifold with 
density (J\f, g, (f>) so that Ric^ = Xg, Asl, therefore we identify Ricci solition with (£, g, <j>, A). 
So, we say that the gradient Ricci solition (£, g, <f>, A) is shrinking, steady or expanding if A > 0, 
A = or A < respectively. 

The condition RicT = Xg says 



therefore, gradient Ricci solitions can be thought as a natural generalization of Einstein metrics. 
One important consequence of F. Morgan result [32] is that gradient shrinking solitions are of 
finite type, which is the key condition to apply Theorem 3.3. 

Lemma 5.1. Let (£,</, 0) be a complete manifold with density so that Ric? 3 > A, A > 0. Then 



In particular, (E,g,<fi) is of finite type. 

Now, we are ready to establish the main theorem in this section. Here, we extend a recent 
result of P. Petersen and W. Wylie [36] relaxing the hypothesis to a L 2 — type bound of the scalar 
curvature R. Moreover, they assume that the Weyl tensor vanishes identically, we change that 
condition by an inequality between the scalar curvature and the Ricci curvature. 



L := A + g(V<f> c ,V-) + \A\ 2 , 



Ric - V (f> = Xg 




dv < +oo. 
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Theorem 5.1. Let (E,g,(f),X) be a complete shrinking soliton of dimension n > 3 so that 
XR < |Ric| 2 . Assume that there exists a sequence of cut-off functions {V'ili ^ Qf^) sucn that 
< ipi < 1 in S, i/te compact sets := /orm an increasing exhaustion of S, and iae 

sequence of weighted energies {J s |V"0i| 2 -R 2 e^ ^s}j is bounded. 
Then, (S, a) aas constant scalar curvature and XR = |Ric| 2 . 

Proof The proof follows from the following equation developed in [36] : 

AR + g(V<p,VR) - 2XR = -2|Ric| 2 . (5.1) 

First, since the scalar curvature satisfies the boundedness condition, R > or (£,<?) is flat 
(see [37]). So, we can continue assuming R > 0. Second, since XR — |Ric| 2 > 0, (5.1) implies 
that 

AR + g(V<f>, VR) = 2XR-2 |Ric| 2 < 0. 
Thus, Theorem 3.3 implies that R is constant, and so, XR = |Ric| 2 . 

□ 

We also have the following classification result: 

Theorem 5.2. • A complete gradient expander soliton of finite type and nonnegative scalar 
curvature is Ricci flat. 

• A complete gradient steady solition of finite type whose scalar curvature does not change 
sign is Ricci flat. Moreover, if (p is not constant then it is a product of a Ricci flat manifold 
with R. Also, S is diffeomorphic to W 1 . 

• A complete gradient shrinking soliton of nonpositive scalar curvature is Ricci flat. More- 
over, (S,g) is isometric to W 1 endowed with the standard Euclidean metric. 

Proof. The proof follows from (5.1). 

• Applying Theorem 3.3 to u = R we obtain that R vanishes identically on S, which implies 
that (£,g) is Ricci flat from (5.1). 

• As above, R satisfies (5.1). So, taking either u = —R if R < or u = R if R > 0, Theorem 
3.3 implies that R = 0. Therefore, (S,g) is Ricci flat. The last sentence follows from [37, 
Proposition 4]. Moreover, S is diffeomorphic to IR n by [6, Proposition 5.7]. 

• Taking u = —R, Theorem 3.3 implies that R = 0. Therefore, (£,g) is Ricci flat. Thus, [33, 
Lemma 7.1] implies that (X,<?) is isometric to M. n endowed with the standard Euclidean 
metric. 

□ 
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6 Weighted stable surfaces and Schrodinger operators 

As we can see from Lemma 1.1, a weighted H^— stable surface has a gradient Schrodinger 
operator, given by (2.1), associated to it. Nevertheless, we can associate a self-adjoint 
Schrodinger operator to a weighted H^— stable surface which is stable in the sense of Definition 
6.1. In this section we will work with surfaces. 

When a Schrodinger operator L := A+q is self-adjoint, we can associated a bilinear quadratic 
form to L given by 

Q : C °°(S) x C °°(S) -> R 

(6.1) 

(u,v) -)• J s (g(Vu,Vv) - quv)dv-z, 

and so 

(-Lu,v) L 2 = Q(u,v) = (-Lv,u) L 2, 
therefore L is self-adjoint with respect to the L 2 — metric. 

Definition 6.1. We say that the Schrodinger operator L is stable if Q(u,u) > for all u G 
C§°(E). 

Now, we will see how to associate a self-adjoint Schrodinger operator to a weighted H$— stable 
surface which is stable in the sense of Definition 6.1. 

Lemma 6.1. Let £ C (N 3 , g , (f) be a complete weighted H^— stable surface. Then, 

[(V- aK)f 2 dv^ < [ \Vf\ 2 dv^ (6.2) 
for any f G Cq°(S) compactly supported piecewise smooth function, where 

v:= \ + \ h * + \\ a \ 2+ \\ v<t) \ 2 ) md a ■= l- ^ 

In other works, the Schrodinger operator 

L:=A-aK + V 

is stable in the sense of Definition 6.1. 

Proof. A tedious but straightforward computation shows that 

Ric^(iV, N) + \A\ 2 = ±R™-K + divpO + F, (6.4) 
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where K is the Gaussian curvature of E, div is the divergence operator on E, and 



1 

2 



F = I(tf 2 + L4| 2 + |v^| 2 ). 



We can also easily check that 
1 



Kg -K + div(X) + F = 3V - K + div(X) + -|V^r. 

Let / e C °°(E) and set u := <p~ 1/2 f £ Cg°(E). Pluging u = .f^ 1 / 2 in (1.3) and using the 
above equalities, we have 

J^V-K + div(X) + ||V0| 2 ^) / 2 ^ s < £ (jVif - / 5 (V/, V0) + ^|V</>| 2 ) cfo s . 

Now, the Divergence Theorem yields 

/ f 2 dw(X)d Vl: = -2 f fg(Vf,V4>) d«s, 

and using the inequality 



fg(Vf,V4>)<2\Vf\ 2 + l\V4>\ 2 f 



we obtain 



j^V-K + ||V0| 2 ^) / 2 ^ s < (V/| 2 + ||V0| 2 / 2 ) ^s, 

or equivalently, 

f (V - aK)f 2 d Vi: < [ |V/| 2 efo s , 

JT, JT, 

as desired. □ 

Recently, the theory of stable Schrodinger operators has been extensively developed since 
the pioneer work of T. Colding and W. Minicozzi [12]. Next, we recall a Colding-Minicozzi type 
inequility (see [7, 21, 29]) we can derive from Lemma 6.1 for a particular choice of radial cut-off 
function. 

As above, we denote by E a connected Riemannian surface, with Riemannian metric g, and 
possibly with boundary <9E. Let po € E be a point of the surface and D(po, s), for s > 0, denote 
the geodesic disk centered at po of radius s. We assume that D(po, s) n <9E = 0. Moreover, let r 
be the radial distance of a point p in D(po, s) to pq. We write D(s) = D(po, s) if no confussion 
occurs. 
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We also denote 

l(s) = Length(d D(s)), 
A(s) = Area(L>(s)), 

K(s) = [ Kdw, 

Jd(s) 

x(s) = Euler characteristic of D(s), 
where length and area are measured with respect to the metric g. 

Lemma 6.2 (Colding-Minicozzi stability inequality). Let £ be a Riemannian surface possibly 
with boundary <9£. Let us fix a point po £ £ and a positive number < s such that D(s)ndT, = 0. 
Assume the Schrodinger operator 

L := A - aK + V 

is stable in the sense of Definition 6.1, where V £ C°°(S) and a is a positive constant, acting 
onfe C °°(£). 

Fix b > 1 and let f : D(s) — > R be the nonnegative radial function given by 

m Ji}-^ * (0.5) 
{ if p£Z\D(s) 

Then, the following holds: 

I Vf(pfdvx < 2anG(s) + mi - b tf + 2a) A(s/2), (6.6) 
JD(s) z s 

where 

G(s) :=- f\f{r) 2 )' X {r)dr<\, 
Jo 

Proof. Since L is stable, plugging / given by (6.5) in Lemma 6.1, we get 

/ Vf 2 dvx< [ (\Vf\ 2 + aKf 2 )dvx. 
JD(s) Jd(s) v ' 

So, from Co-Area Formula and Fubbini's Theorem, we obtain 



/ Vf 2 dv s < [ (\Vf\ 2 + aKf 2 )dvj; 
JD(s) JD(s) v ' 

= [ \f{r) 2 I ldl + af{r) 2 [ K{r)dl\dr 

JO \ JdD(r) JdD(r) J 

= [ {f'{r) 2 l(r) + af(r) 2 K'(r))dr 
Jo 
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Now, from [42] and since (f(r) 2 )' = 2/(r)/'(r) < 0, we have 

-(f(rf)'K(r) < (/(r) 2 )'(/'(r) - 27T X (r)). 
Integrating by parts and taking into account that J Q s (f(r) 2 )' = — 1, we obtain 

/ Vf 2 dw < -2na f (f(r) 2 )' X {r)dr + f (a(/(r) 2 )'/'(r)) dr 
Jd(s) Jo Jo 

= 2irG(s) + f ((1 - 2a) f (r) 2 - 2af(r)f"(r)) l(r)dr 
Jo 



26-2 

l{r)dr 



< 2lG{s) + "(Mi-y + 2» )j4 ° V2) 

Note that the bound on G(s) follows since the Euler characteristic of D(s) is less than or 
equal to 1. □ 

Before we finish this secion, we will see an useful Lemma: 

Lemma 6.3. Under the conditions of Lemma 6.2, i/S is complete and there exists sq such that 
X(s) < -M, M > 0, for all s > s , then 

G(s) < -(M + l)/( So ) 2 + l. 

Proof. Assume there exists sq so that for all s > sq, we have x( s ) ^ —M. Therefore, following 
Lemma 6.2, we have 

G(s) = - / (f(r) 2 Yx(r) = - / (f(r) 2 )'x(r) - / (f(r) 2 )'x(r) 
Jo Jo Js 

< - / (f(r) 2 )' + M / (f(r) 2 )' = - (f(s ) 2 - f(0) 2 ) + M (f(s) 2 - f(s ) 2 ) 

JO Js 

= -(M + l)/( So ) 2 + l, 

since —{f{r) 2 )' > and x( r ) < 1 f° r an r - D 

7 Manifolds with nonnegative Perelman Scalar Curvature 

In this Section, we will study complete weighted H^— stable surfaces in three-manifolds with 
density (Af,g,4>) and nonnegative P-scalar curvature > 0. 

The first thing we shall observe is when we can have a complete noncompact weighted 
H^— stable surface under conditions on the P-scalar curvature. The following result is the 
extension of the diameter estimate for stable H— surfaces given by H. Rosenberg 
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Theorem 7.1. Let S C (Af,g,(f>) be a weighted H^— stable surface with boundary d'E. If K£ + 
#| > c > on S, i/ten 

27T 

cfe (p, 9S) < -7= /or a// p G S, 
V3c 



where d^, denotes the intrinsic distance in £. Moreover, if £ is complete without boundary, then 
it must be topologically a sphere. 

Proof. The condition i?^ 3 + > c > implies that V > c > for some positive constant c. 
Therefore, [29, Theorem 2.8] with a = 1/3 implies the result. We will sketch the proof for the 
reader convinience. 

Since L is stable, there exists a positive function u such that Lu = [22]. Consider the 
conformal metric g := u G g, where g is the metric on S. Denote by K and X the Gaussian 
curvature of g and g respectively. So, since = 0, we get 

*_ 3.^ (v + !^)>S^(c (7.1) 

Take p G S and let 7 be a 5— geodesic ray emanating from p. Denote by Z and I the length of 
7 with respect to 5 and g respectively. Since 7 is a g— minimizing geodesic, the Second Variation 
Formula of the arc-length gives that 



i ( 



) ds > 0, (7.2) 



for any smooth function : [0,Z~] ->• R such that 0(0) = 0(f). From (7.1), (7.2), |Vu| > 
(it o 7)'(s) = u'(s) and changing variables ds = u 3 ds, we get 

f (V (*) 2 - 3 [c + ^) 0(s) 2 ) ds > 0. (7.3) 

Take (f> = u 3 / 2 ^, where ip : [0, 1] — > M. is a smooth function such that ip(0) = ip(l) = 0. Then, 
the above (7.3) yields 

f (-4i}"(s)*p(s)-3cip(s) 2 )ds>0, (7.4) 
Jo 

for all V G Cg°([0,i]). Taking ip(s) = sin (fs) in (7.4), we get 



l / Ji 







4^r — ;>r I sin 



ys) ds > 0, 
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which implies 

l < *J± , (7.5) 

which gives the desired estimate. 

Now, if £ is complete, then (7.5) and the Hopf-Rinow Theorem imply that £ must be 
compact. Moreover, applying the operator L to the test function 1, we have 

\ [ KdvT,> cArea(E), 

which implies, by the Gauss-Bonnet Theorem, that x(£) > 0. 

□ 

The above result says that, in a manifold with density (Af, g, 4>) and K£ > 0, the only 
complete noncompact weighted -ff^— stable surface we shall consider are the 0— minimal ones. 
So, the next step is to extend the well-know result of D. Fischer-Colbrie and R. Schoen [23] on 
the topology and conformal type of complete noncompact stable minimal surfaces. P. T. Ho [17] 
extended Fischer-Colbrie-Schoen result under the additional hypothesis that |V</>| is bounded 
on £. 

Here, we drop that conditon and caracterize the weighted stable cylinder as totally geodesic, 
flat and the density must vanish along £: 

Theorem 7.2. Let S C (N 3 ,g,(j)) be a complete (noncompact) weighted stable minimal surface 
where > 0. Then, E is conformally equivalent either to the complex plane C or to the 

cylinder E> 2 x R. In the latter case, E is totally geodesic, flat and 4> is constant and = 
along E. 

Proof. From Lemma 6.1, L := A — aK + V is stable in the sense of operators with 1/4 < a = 1/3 
and V > 0. From [7] or [29], we can see that E is conformally equivalent either to the plane or 
to the cylinder. 

In the latter case, either [19, Lemma 2.1] or [3, Theorem 1.3] give us that K = and V = 0, 
which implies that \A\ 2 = and <fi is constant and = along E. 

We will sketch the proof for the reader convinience. Let us consider the radial function 

J(l - r/sf r<s 
/(r):= | r> S ' 

where r denotes the radial distance from a point po £ S. Then, from Lemma 6.2, we have 

/ (1 - r/s) w V dw < lirG(s) - \a(s/2). (7.6) 

JD(s) 3 S 
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First, since E is complete, from (7.6) and Lemma 6.3, we get that x(E) equals either to or 
1, that is, £ is topologically either the plane or the cylinder. 
Second, again from (7.6) and V > 0, we obtain 

for some constant C, which yields that E has quadratic area growth. Thus, since £ has finite 
topology and quadratic area growth, each end of £ is parabolic (see [12]). Therefore £ is 
conformally equivalent either to the plane or to the cylinder. 

Now, suppose that E is conformally equivalent to the cylinder. We will show that E is flat, 
totally geodesic and <p is constant and K£ = along E. 

• Claim A: V vanishes identically on E. That is, E is totally geodesic and 4> is constant 
and = along E. 

Suppose there exists a point po £ Y> so that V(po) > 0. From now on, we fix the point 
Pq. Then, there exists e > so that V{q) > 5 for all q £ D(e) = D(po,e). Since E is 
topollogically a cylinder, there exists so > so that for all s > so we have x( s ) ^ (see 
[7, Lemma 1.4]). 

Now, from the above considerations and (7.6), there exists (3 > so that 

< p < 2ttG(s). 

But, from Lemma 6.3, we can see that 

G( s ) = -/(s ) 2 + l = -(l- V*) 2 + l, 

therefore, 

G(s) < 1 — (1 - s /s) 2 -)• 0, as s ->■ +oo, 
which is a contradiction. Thus, V vanishes identically along E. 

• Claim B: E is flat. 

Since E is a cylinder and (7.6), we have 

lim ^1=0, 

therefore, K > on E (see [3, Proposition 5.1]). Thus, Cohn-Vossen Theorem yields that 
K vanishes identically on E. 

□ 
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7.1 Stable hypersurfaces 



We finish this section classifying stable hypersurfaces in manifolds with density that verify a 
lower bound on its Bakry-Emery-Ricci tensor. Specifically: 

Theorem 7.3. Let £ C (ftf, g, (f>) be a complete H^— stable hypersurface of finite type and assume 
that Ric^ 1 > k, k > 0. Then, £ is totally geodesic and Ric^(iV, N) = along £. 

Proof. Since £ C (J\f,g,(j>) is .ff^— stable, the gradient Schrodinger operator, given by (2.1) 

L^u := Au + g (V<p, Vu) + q^u, 

where q^ := Ric^ D (A r , N) + |A| 2 , is stable in the sense of Definition 6.1. Therefore, from Lemma 
3.1, there exists a positive subsolution, i.e., there exists u > such that L^u < 0. Moreover, 
q^ > since Ric^ > 0. 

Thus, Theorem 3.3 implies that q^ = 0, that is, £ is totally geodesic and Ric^°(iV, N) = 
along £. □ 



8 Local splitting 

In this section we study the topology and geometry of a manifold with density (J\f, g, 4>) verifying 
some lower bound on its P-scalar curvature R£ ■ The idea is to extend the splitting theorems 
developed by Cai-Galloway [5], Bray-Brendle-Neves [4] and I. Nunes [34]. Here, we will take the 
unified point of view considered by Micallef-Moraru [31]. 

First, we will generalize an area estimate for compact stable minimal surfaces in three- 
dimensional manifolds with a lower bound on its scalar curvature to the context of weighted 
stable minimal surfaces: 

Proposition 8.1. Let £ C (M,g,4>) be a compact weighted stable minimal surface. 

1. If R£ > Xe^, for some positive constant X, then A^,(£) < Moreover, if A^iTi) = 

then £ is totally geodesic, the normal Barky- Emery- Ricci curvature of (jV, g, 4>) vanishes 
and cf> is constant along £. Moreover, R^ = 2K = Xe^ along £. 

2. If R^ > and £ has genus one, then £ is totally geodesic, the normal Barky- Emery- Ricci 
curvature of (M,g, 4>) vanishes and <f> is constant along £. Moreover, R^ = 2K = along 
£. 

3. If R^ > —Xe^, for some positive constant X, and^f = genus(£) > 2, then A^(T,) > 47r( ^ 1 ) , 

Moreover, if A^T,) = 47r (^ -1 ) > then £ is totally geodesic, the normal Barky-Emery-Ricci 
curvature of (N,g : 4>) and vanishes and (p is constant along £. Moreover, R^ = 2K = 
—Xe^ along £. 
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Proof. Since E is compact and weighted stable, from Lemma 6.1 taking / = 1 in (6.2), we get 

which yields, from Gauss-Bonnet formula, the following 

\ [ ^d« E <47r(l- 7 ), (8.2) 

where 7 := genus (E). 

1. If > Ae^, for some constant A > 0. 
From (8.2), we have 



^(S) = ±Je*dvz<ljR$' dvx < 4tt, 



note that E must be a topological sphere. Therefore, from this last equation, the inequality 
follows. 

Moreover, in the case equality holds, (8.1) implies that E is totally geodesic and <p is 
constant along E. This implies that the Schrodinger operator associated to E given by 
Lemma 6.1 reads as 

L:=A + ±(R^-2K). 
Therefore, the equality implies that 

{Rf - 2K)dvx = 0, 



which yields that the first eigenvalue of L is zero, i.e., Ai(— L) = 0, and so, the constant 
functions are in the kernel of L, therefore 

Rf = 2K along E. 

Thus, from (6.4) and the fact that \A\ 2 = and (j> is constant along E, we obtain 

Ric^(iV, N) = Ric^(iV, N) + |^| 2 = ±fl£> - K = along E, 
which finishes the proof of item 1. 
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2. If > and E has genus one. 

As above, taking / = 1 in (8.1), we get 

which implies that E is totally geodesic, (p is constant along E and = along E. So, 
arguing as above, we get K = and Ric^ 1 (N, N) = along E as well. 

3. If > —Xe^, for some constant A > 0. 
This case is completely analogous to item 1. 

□ 

Next, we shall prove the existence of a one parameter family of constant weighted mean cur- 
vature surfaces in a neighborhood of a totally geodesic compact surface verifying the conditions 
on Proposition 8.1. 

Proposition 8.2. Let E C (N 3 ,g,(fi) be a compact immersed surface with unit normal vector 
field N. Assume that 

• E is totally geodesic and (ft— minimal, 

• the normal Barky-Emery-Ricci curvature of (N,g,(j)) vanishes and <p is constant along E. 

Then, there exists e > and a smooth function w : E x (—e,e) — > R such that, for all 
t £ (—e,e), the surfaces 

E t := {e Wp (w(p,t)N(p)) : p £ E} 
have constant weighted mean curvature Hj,(t). Moreover, we have 



w(p,0) = 0, ^ 



w(p, t) = 1 and I (w(-, t) — t)dv-£ = 0, 
=o Jx 

for all p £ E and t £ (— e, e) . 

Proof. We follow [34]. Fix a £ (0, 1) and consider the Banach spaces 

X[n) := ju G C n ' Q (E) : J ud Vj] = oj for each n £ N. 

For each u £ C 2 ' a (E) we define 

S(u) := {exp p (n(p)iV(p)) : p £ E} . 
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Choose £\ > and 5 > so that E(u + i) is a compact surface of class C 2,a for all (t,u) G 
(-ei,ei) x 5(0,5), here 5(0,5) := {u G C 2 ' Q (E) : |«| c2 ,« < 5}. Denote by 5^(u + t) the 
weighted mean curvature of E(u + i). 

Consider the map $ : (-£i,£i) x (5(0, 5) n X{2)) X(0) given by 



$(t,u) :=ffy(u + t)- 



5 (u + i)e^ cfo E . 



First, note that $(0,0) = since E(0) = E. Second, we compute 5<I>(0,0) • v for any 
v G X{2). We have 



D$(0,0) -u 



ds 



$(0,sw) = 



s=0 



(is 



1 



rf **<*"> " .4,(E) h i> 



d 



H^sv) dvx 



L * v - am L {L * v) e " =Av ~a^)L Av 



s=0 
dVT, 



= Av 

since 4> is constant along E and 



d_ 

ds 



H^sv) = L^v = Av 



s=0 



from (2.1) and the hypothesis. 

So, since A : X (2) — > X(0) is a linear isomorphism, by the Implicit Function Theorem, there 
exists < e < £i and u(t) := u(t, •) G 5(0,5) for t G (—£,£) such that 



u(0) = and $>(t,u(t)) = for all t G (-£,£). 



Finally, defining 



w(t,p) = u(t,p)+t, (t,p) G (-£,£) X E, 

we obtain the result. 

Now, we are ready to prove the local splitting result: 



□ 



Theorem 8.1. Let (N 3 , g , (f) be a complete manifold with density containing a compact, em- 
bedded weighted area-minimizing surface E. 

1. Suppose that > Xe^, for some positive constant X, and ^(E) = Then E has genus 
zero and it has a neighborhood in M which is isometric to the product § 2 x (— e, e) with the 
product metric g + \ + dt 2 (up to scaling the metric g); here g + \ is the metric of constant 
Guassian curvature +1. Moreover, <f> is constant inU. 
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2. Suppose that > and E has genus one. Then E has a neighborhood in M which is flat 
and isometric to the product T 2 x (— e,e) with the product metric go + dt 2 ; here go is the 
metric of constant Guassian curvature 0. Moreover, <j) is constant inU. 



3. Suppose that R^ > — Ae^, for some positive constant X, E has genus 7 > 2 and A^iT,) = 
47r ^ — . Then E has a neighborhood in M which is isometric to the product E x (— e,e) 
with the product metric g-\ + dt 2 (up to scaling the metric g); here g_\ is the metric of 
constant Guassian curvature — 1. Moreover, (p is constant inU. 

Proof. Let E C (J\f, g, 4>) be a compact embedded weighted area-minimizing surface under any 
of the conditions above. From Proposition 8.1 and Proposition 8.2, there exists e > and a 
smooth function w : E x (— e,e) — > K. such that, for all t G (— e, e), the surfaces 

£* := {exp p (w(p,t)N(p)) : p G E} 

have constant weighted mean curvature H^it). Moreover, we have 

d 



w(p,0) = 0, ^ 



w(p, t) = 1 and / (w(-,t) — t)dv-£ = 0, 
=0 



for all p G E and t G (— e, e). 

On the one hand, from (2.1) we have: 



L^u = Au + (Ric^(iV, N) + \A\ 2 )u + g (V0, V«) 

= Au + (^R$> - K + divX + F)u + g Vn) 

>A.+ (i^-K + divX + J F)n-^ |Vu| 



that is, 



Let 



2 9 ' 2 2n 

A.-L^ + ^-K + divX + F-L^). 

nAlnn+ + (1^ - K + divX + 5 + J4!) u> 



IV7/I 2 1 #i Ml 2 

L n > uA In u + L_L + (-R°° - K + divX + ^ + U-) u . (8.3) 



/t(x) = exp p (w(p,t)N(p)), p G E, i G (-e,e), 
thus /t(E) = E t for each t G (— £, e) being E = S. Define the lapse function p t : E — > R by 

ft(p)=5(iVt(p),^/t(p) 
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where N t is a unit normal vector along Y, t . Set £ : £ x (— e, e) — > Af given by 

= ft(p), (P,t) G S X (-£,£■). 

By the First Variation Formula for the Generelized Mean Curvature (see [2] and [8]), we 
have 

Hfo) = -L^pt. 
On the other hand, from (8.3), we obtain 

H'^t)- < -A t In ^ - d\wX t + K t - \r$. (8.4) 
Pt 2 

1. If i?oo > Ae^, for some positive constant A, and A,p(E) = 

Since i?^ 1 > Ae^, integrating (8.4) over £ t , we obtain for any i G (0,e) 

/ -d« St = / ^)-^s t < / ^^s t -^ / e^ S( =47r-^(i), 



that is, 



H'4t) f }-dv^<±(A^)-A^t)) 
J Ht Pt * 



where we have used that H^it) is constant on each S t . Since po = 1 on S, by continuity, 
there exists e > so that 1/2 < p t < 2. This yields that 



JUt Pi 1 



since S is a minimizer for the weighted area funtional. Therefore, H'^t) < for all 
£ G (0, e). So, by the First Variation Formula for the weighted area, we have 

A'^t) = [ H^{t)p t e* dv^ t < for all t G (0, e), 

and using that S is a miminizer for the weighted area, we obtain A^it) = A^iO) for all 
t G (0, e). Analogously, we can prove that A^it) = A^fO) for all t G (— e, 0). 

Thus, Proposition 8.1 implies that is weighted minimal, totally geodesic, <j) = <f>t is 
constant along S t and -R^ 1 = 2K t = Ae^* (constant) along £ for each t G (— e, e). Now, 
since S t is weighted minimal and totally geodesic for each i G (— e,e), we get 

<7(Vln^AT t ) =0, 

which yields that (j) = (fro is constant in U := £(X x (— e, e)) C M . Moreover, from (8.4) 
and St being totally geodesic, we obtain that pt is constant along each Sj, t G (— e,e), i.e., 
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pt is a function of t only. This yields that the vector field N t is parallel on U which implies 
that the integral curves of N t are geodesies and so 

£(t,p) =exp p (tN(p)) for all p € E, 

futhermore, the above map is an isometry from E x (— e, e) to ZY, where we consider the 
product metric + dt 2 in E x (— e,e). So, scaling g, we can assume that K t = 1 for all 
E t , which finishes item 1. 

2. If iioo > and E has genus one. 

The proof is completely analogous to the one on item 1. 

3. If i?oo > — Ae^, for some positive constant A, E has genus 7 > 2 and A^(E) = inili. 
As we did in item 1, integrating (8.4) over E t , we obtain for any t £ (0, e) 

[ H'^dv^K I K t dw t + \l e*d^ t = 27r(l- 7 ) + ^(t), 



that is, 



flj(t) / 1 dtfc, < £ (^(t) - A^O)) = ±f A'4s) ds, 
J Ht Pt z * Jo 



so, using the first variation formula for the weighted area and that H${t) is constant on 
each Et we get 

H'+it) I -dw t < f H^{s)( I p s dvAds. (8.5) 

JT, t Pt JO \Jt. s J 

Assume there exists to £ (0,e) so that H^to) > 0. Define 

I:={te [0,to] : H^t) > H<p(t )} . 

We will show that inf / = 0. Assume by contradiction that inf I = i > 0. First, rewrite 
(8.5) as 



where 



H'^K-^J H^s^ds, (8.6) 

ip(t) := / — dv£ t and fx{t) = [ pte^ dvs t . 
Jt,i Pt JT, t 
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Since po = 1 and <f> = 4>q on S, by continuity, there exists £ > so that 1/2 < p t < 2 and 
< < 2e^°. This yields that 

< M (t) < 2A (f> (t) , for all t G (-e,e) 

and 

^ < < 4e*>^(t) , for all t G (- e ,e) 

Shrinking e if necessary, we can assume A ( f,(0)/2 < < 2A<f,(0) for all t G (— e, e). 

Therefore, we obtain 



< and //(i) < 4^(0) for all t G (-e,e) 



^(*) ^(0) 

Second, by the Mean Value Theorem, there exists t' G (0, t) so that 

H^t) = H'^t')i. (8.7) 

Therefore, from (8.6) and (8.7), we get 

H^t') <^fjf Q H 4> {s)p{s)ds < 32e*°ffy(t / )t'* < S2e^e 2 H (j> (t'), 
so, if e < , 1 , , we get the desired contradiction. 

Since inf / = 0, it follows that H^O) > H^it) > 0, which contradicts that £ is weighted 
minimal. Therefore, H$(t) < for all t G (0, e). Now, we can do the same for t G (— e, 0). 
And we can finish as in item 1. 



□ 



9 Topology of three-manifolds with density 

This section is mostly based on Section 8 and the following observation: 

Proposition 9.1. Let £ C (ftf 3 ,g,(fi) be a compact surface in a manifold with density. Then, 
S is weighted area-minimizing if and only ifH is area-minimizing in (J\f,g := e^g). 

First, we extend the rigidity result of Bray-Brendle-Neves [4]: 
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Theorem 9.1. Let (N 3 ,g,(f)) be a compact manifold with density so that is positive and 
tt 2 (M) / 0. Define 

A(AT,g,4>) :=inf{A (/(§ 2 )) : f € ?} , 

where T is the set of all smooth maps f : S> 2 — > N which represent nontrivial elements in TT2(M). 
Then, we have 

A(N, g, 0)min } : x G N \< 8tt. (9.1) 

Moreover, if equality holds, the universal cover of (M,g) is isometric (up to scaling) to the 
standard cylinder § 2 x R and <f> is constant in M . 

Proof. Since is positive and N is compact, we get 

f R°°(x) ] 
> Ae^, where A := min <^ ^ : x G N \ . 

So, from Theorem 8.1, we get the upper bound. 

Now, assume the equality holds in (9.1). From Propostion 9.1 and results of Meeks-Yau [30], 
there exists a smooth inmersion / G T so that 

A <i> {fi^))=A{N,g,4>). 

Denote X = /(§ 2 ). Since we are assuming equality, Theorem 8.1 asserts that locally, in a 
neighborhood of S, Af splits as a product manifold and <p is constant. We can also see that each 
leaf of the product structure is weighted area minimizer, so we can continue the process and we 
construct a foliation {Xt} tg ]g>, 



E t := {exp /(p) (tiV(p)) : p G s| 



of embedded spheres which are totally geodesic, has the same constant Gaussian curvature and 
4> and R^ are constant. 

Now, it is not hard to see that the map 

$ : § 2 x R -»• Af, 

given by 

$(p, i) = exp /(p) (tiV(p)) 

is a local isometry (see [4, Proposition 11]). From here, it follows that $ is a covering map and 
therefore, the universal cover of (M, g) is isometric to S 2 x M equiped with the standard metric. 
This finishes the proof of Theorem (9.1) □ 
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Second, the rigidity result of Cai-Galloway [5]. The proof is completely similar to Theorem 
9.1 so, we omit the proof. 

Theorem 9.2. Let (N 3 ,g, (ft) be a complete manifold with density so that is nonnegative. 
If (M, g, (ft) contains a weighted area miminizing compact surface in its homotopy class of genus 
greater than or equal to 1, then the product manifold T 2 xM, where T 2 is a torus equiped with the 
standard flat metric, is an isometric covering of {M ,g), and (ft is constant in M . In particular, 
{M, g) is flat. 

And finally, the rigidity result of Nunes [34] : 

Theorem 9.3. Let (AA 3 , <?, (ft) be a complete manifold with density so that R™ > — Ae^ for some 
positive constant A. Moreover, suppose that £ C (M,g,4>) is a two-sided compact embedded 
Riemannian surface of genus 7 > 2 which minimizes area in its homotopy class. Then, 

ME) > M^i). 

Moreover, if equality holds, the universal cover of (N,g) is isometric (up to scaling) to the 
product manifold S 7 x R, where £ 7 is a compact surface of genus 7 equiped with a metric of 
constant Gaussian curvature —1, and (ft is constant in N . 



10 Further comments: Optimal Transportation 

As we said at the Introduction, given a manifold with density (M,g,4>), (ft £ C°°(M), one can 
consider eft as a function to distorted the reference volume meassure dv, i.e., dm = dv. So, we 
can consider the Boltzman entropy functional 

H ra (fi) := / p\n pdm , /i = p = pdm £ T>2(N), 

where ^(AA) is the Wasser stein space. 

Following Otto's calculus, we saw that critical points of the Boltzman entropy, i.e., p G 7^2 (AA) 
so that 

grad M "H m = 0, 

correspond to positive solutions to the gradient Schrodinger operator 

L^p:= Ap + g(V(ft,Vp). 

Now, since optimal transportation plans (for the Boltzman entropy functional) must be 
critical points, we call them critical transportation plans, we can locate the space of optimal 
transportation plans under global conditions on the initial data (N,g,(ft) using Theorem 3.2, 
that is: 
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Let (Af,g,(f>) be a complete manifolds with density of finite type. Then, critical trans- 
portation plans fj, for the Boltzman entropy functional are := a dm, a G M + . 
Moreover, if we assume 

/ dfi = dm, 

then, the only optimal transportation plan is fi := 4>dv. 

Actually, this is the situation for the Gaussian meassure space (R n+1 , (•,•), e -71 "^ 2 ). Note 
that, if dm := e^ dv € V2(N), then (J\f, g, 4>) is of finite type. 

Remark 10.1. It would be interesting to investigate deeper this relationship, even for other 
functional and not only for the Boltzman entropy. 
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